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ly-^ ' Abstract 
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' The metric for plane gravitational waves is quantized within the 

^ ■ Hamiltonian framework, using a Dirac constraint quantization and 

G^ . the self-dual field variables proposed by Ashtekar. The z axis (direc- 

I tion of travel of the waves) is taken to be the entire real line rather 

CT' than the torus (manifold coordinatized by (z,t) is RxR rather than 

Si X R). Solutions to the constraints proposed in a previous paper 
involve open-ended flux lines running along the entire z axis, rather 
I than closed loops of flux; consequently, these solutions are annihi- 

^ ■ lated by the Gauss constraint at interior points of the z axis, but not 

■ at the two boundary points. The solutions studied in the present 

paper are based on closed flux loops and satisfy the Gauss constraint 
for all z. 
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1 Introduction 

This is one of a series of papers discussing quantization of plane grav- 
itational waves [|T], 0. The emphasis is on the Ashtekar approach 
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0, |]; a major goal is to understand how radiative phenomena are 
expressed using the Ashtekar language. The method of quantization 
is canonical. Due to the high degree of symmetry (full translational 
symmetry in the x,y directions), four of the seven scalar, vector, and 
Gauss constraints can be solved and eliminated from the problem 
before quantization. (The four are: the two vector constraints gen- 
erating x,y spatial diffeomorphisms; and the two Gauss constraints 
generating "internal" rotations around the X and Y axes.) The re- 
maining three constraints are imposed after quantization, and the 
key step is finding a wavefunctional which is annihilated by these 
constraints. A previous paper (II) wrote down the threee survivng 
constraints in the Ashtekar formalism, quantized them by replacing 
one half of each (A, E) canonical pair by a functional derivative, 
and wrote down a class of wavefunctionals annihilated by the the 
constraints 0. 

The ip constructed in II are unusual. It is usual to have the ip 
depend on the connection fields A, in order to exploit analogies with 
Yang-Mills theory, whereas these depend on only one connection 
field, A^ • The remaining dependence is on densitized inverse triad 
fields E. Each ip is a string of E operators separated by holonomy 
matrices depending on A^. If the holonomies are visualized in the 
usual way as lines of flux, tjien these solutions are flux lines extend- 
ing along the z axis, with E fields inserted into the line. 

These solutions satisfy the Gauss constraint at finite z, but not 
at infinity. This can also be understood in visual terms. The flux 
line is an open, rather than closed curve, since the flux line does not 
loop back to close on itself. Since flux exits at infinity, the Gauss 
constraint is not satisfied there. 

The present paper constructs solutions which are "closed flux" 
in nature and are annihilated by the Gauss constraint everywhere. 
Since the present paper is for the most part an extension of II, I have 
refrained from repeating many matters of philosophy and motivation 
already discussed in that paper. I have repeated here some of the 
more detailed results of II, however, so as to make the calculations 
of the present paper self-contained: see the first half of section 2 and 
the opening paragraphs of section 3. 

The first half of section 2 reviews the quantization procedure of 
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paper II and writes out the three constraints using Ashtekar vari- 
ables. The second half of section 2 constructs a point-splitting regu- 
larization for one of the terms in the scalar constraint. (In II it was 
not necessary to regulate this term, because it vanishes identically 
when acting on the solutions of II.) Section 3 presents the ansatz 
for and shows that the three constraints annihilate the ansatz. 
Section 4 presents conclusions and directions for further research. 

My notation is typical of papers based upon the Hamiltonian 
approach with concomitant 3+1 splitup. Upper case indices A, B, 
. . .,1, J, K, . . . denote local Lorentz indices ("internal" SU(2) indices) 
ranging over X, Y, Z only. Lower case indices a, b, . . ., i, j, . . . 
are also three- dimensional and denote global coordinates on the 
three-manifold. Occasionally the formula will contain a field with a 
superscript (4), in which case the local Lorentz indices range over 
X, Y, Z, T and the global indices are similarly four-dimensional; or 
a (2), in which case the local indices range over X, Y (and global 
indices over x. y) only. The (2) and (4) are also used in conjunction 
with determinants; c. g., g is the usual 3x3 spatial determinant, 
while '•^•fe denotes the determinant of the 2x2 X, Y subblock of the 
densitized inverse triad matrix E^. I use Levi- Civita symbols of 
various dimensions: exxYZ = exYZ = cxy — +1- The Minkowski 
metric convention is t/tt = — 1- The basic variables of the Ashtekar 
approach are an inverse densitized triad and a complex SU(2) 
connection A^. 



e^ is the usual 3x3 triad matrix and e^ is its inverse. The local 
Lorentz indices are vector rather than spinor, and strictly speaking 
the internal symmetry is 0(3) rather than SU(2). After the gauge- 
fixing allowed by the planar symmetry, the internal symmetry is 
0(2) rather than U(l). 

2 Quantization 

The first part of this section summarizes the main results and for- 
mulas on quantization of plane waves from II, for the convenience of 




Vgel; 



(1) 




3 



readers who do not have a copy of II at hand. The second part dis- 
cusses a regularization needed to define unambiguously the action 
of one term in the scalar constraint. 

I follow the quantization procedure used by Husain and Smolin 
P]. These authors use the planar symmetry to solve and eliminate 
four constraints (the x and y vector constraint and the X and Y 
Gauss constraint) and correspondingly eliminate four pairs of (E^, 
A^) components. The 3x3 matrix then assumes a block diagonal 
form, with one 1x1 subblock occupied by E| plus one 2x2 subblock 
which contains all E^ with a = x,y and A = X,Y. The 3x3 matrix 
of connections A^ assumes a similar block diagonal form. None of 
the surviving fields depends on x or y. 

After these simplifications, the total Hamiltonian reduces to a 
linear combination of the three surviving constraints, 

Ht = N'p)EeABeab(E|)''A^AB/2 + eMNE^F^,] 

-iNG[a,E|-eijEfAJ] + S.T. 
= N'Hs + N^H, + NgHg + S.T., (3) 

where 

F^b = d.A^ - ^nqA^^AQ (4) 

Hs,Hz, and Hq are the surviving scalar, vector, and Gauss con- 
straints. Strictly speaking these are Hamiltonian densities; for sim- 
plicity I have suppressed an integration over the z axis. S.T. denotes 
"surface" terms (terms evaluated at the two endpoints on the z axis). 
The detailed form of these terms is worked out in II but will not be 
needed here. 

The theory may be quantized by replacing half of each conjugate 
momentum-coordinate pair by a functional derivative, in the usual 
way, 

A^ ^ h5/5E%; A = X,Y;a = x,y; (5) 
E| ^ -h6/6Al (6) 

in order to realize the canonical commutation relations, equation (0) 
In II I modified the Husain-Smolin quantization procedure in 
certain respects. Those authors consider the closed topology (the 



4 



z axis is topologically a circle Si) whereas I consider the open one 
(the z axis is the real line R). Persons familiar with open topology 
calculations in three spatial dimensions might suppose that there- 
fore the constraint equation should be integrated over the entire 
real axis, z = — oo to z = +00. In the present case (effectively one 
space dimension rather than three, because of the planar symme- 
try) the space does not become flat at z goes to infinity. Moving 
the boundary points to spatial infinity does not especially simplify 
matters, and in fact integrals are better behaved if all integrations 
are taken from zq to Z2n+i (in the notation to be used in section 3), 
where the points zq and Z2n+i are a finite distance from the origin. 
The radiation is assumed to be wave packets confined to the region 
Zq < z < Z2n+i iusidc the boundaries, so that the gravitational field 
is relatively simple (but not fiat) at boundaries. Again, the details 
of the fields at the boundary will not be needed here, but perhaps 
it is helpful to point out that the result that the space does not 
become fiat as z goes to infinity does agree with one's intuition from 
Newtonian gravity, where the potential in one spatial dimension due 
to a bounded source grows as z at large z. 

I have modified the Husain-Smolin quantization in another re- 
spect. In equation (^ the Lagrange multiplier N' and the scalar 
constraint Hs are rescaled versions of the usual Lagrange multiplier 
N and Ashtekar scalar constraint H: 

NH = [(NE|)][H/E|] 

= [(N')][Hs] (7) 

A rescaling of this type, for a one space-dimensional theory, was 
first proposed by Teitelboim in the context of geometrodynamics . 
This innocuous-looking rescaling has profound effects on the closure 
of the algebra of constraints. As shown in II, when the factors are 
ordered as in equation (^, with all functional derivatives to the 
right, then the constraint algebra is consistent, and simultaneously 
Hz can be interpreted as the generator of diffeomorphisms. 

However, since Hs = H/E^, and H is polynomial, the new scalar 
constraint Hs is rational. I view this complication 

as a price which must be paid; but it is a small price, consid- 
ering what one gets in return. The constraints close, which they 
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should as a matter of principle; and the quantum theory possesses 
the diffeomorphism invariance characteristic of the classical theory 



Nevertheless, at some point a price must be paid: how does one 
define the inverse operator (E|)~^? Fortunately, for the wavefunc- 
tionals considered in this paper and in II there is a natural defini- 
tion of the inverse. The variable conjugate to E| occurs only in 
holonomies 

M(zi+i,Zi) = exp[i r''' S,A,2(z')dz']. (8) 

Sz is the usual 2j+l dimensional matrix representation of the gener- 
ator for SU(2) rotations around the Z axis, and there is an explicit 
factor i because is Hermitean. From equation the action of E| 
on this function is especially simple, since it merely brings down a 
factor of — iSz = — im^ where m^ is the eigenvalue of S^. The natural 
definition of the inverse (E|)-i is then 

(E|)-iM= (-im,)-^M. (9) 

This works provided m^ is never allowed to assume the value zero. 

So far I have been merely reviewing matters already discussed 
in II. I now consider a new topic, regularization of the constraints. 
As with the definition of the (E|)~^ operator just given, the reg- 
ularization will make the constraints finite and unambiguous when 
the constraints act on the wavefunctionals considered in this paper; 
no claim is made that the proposed regularization will work in all 
circumstances. I follow Husain and Smolin P, who suggest a sim- 
ple point splitting regularization. This simple technique works be- 
cause the planar case is simpler than the full three-dimensional case. 
In the full case, the wavefunctionals typically involve integrals over 
loops (Wilson loops) which are one-dimensional, while the functional 
derivatives analogous to equation give rise to three-dimensional 
delta functions |TD[. There are more delta functions than in- 



tegrals, which is a recipe for divergences. In the planar case, the 
wavefunctionals still involve one-dimensional integrals ( always over 
the z axis), but the delta functions from the functional derivatives 
are also one- dimensional. 
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Even when the dimensionahty of the integrals matches the di- 
mensionahty of the delta functions, it is possible to get a divergence 
if the wavefunctional contains a product of two or more fields evalu- 
ated at the same point. (For example a wavefunctional of the form 
/ dz'[E(z_')]^ acted upon by a constraint operator containing a prod- 
uct [(5/5E(z)]^ gives rise to an undefined [(5(z — z')]^ .) But this does 
not happen in the present case: the wavefunctionals studied here 
involve at most one field evaluated at each point. 

In fact only the first term in equation (Q) needs regularization, 
and the regularization is needed to remove an ambiguity rather than 
a divergence. This term is the only one containing three, rather than 
one functional derivative. The wavefunctional is path-ordered, so 
that one can talk meaningfully about "adjacent" fields. When the 
first term acts on three adjacent fields, the point splitting is needed 
to define unambiguously which three fields will be acted upon in 
which order. 

I call the first term He (E for the *^^-*E which it contains) and point 
split it as follows: 

He = N'(2)E(z)eABeabA^(z + e)[E|(z)]-iA^(z - e)/2. (10) 

One can now see how the point splitting removes an ambiguity. 
It defines unambiguously which holonomy M the [E|]~^ operator 
acts on, wheri the wavefunctional contains a chain of holonomies 
separated by E fields. 

J dzidz2 ■ ■ ■ M(zi+i, Zi)E(zi)M(zi, Zi_i)E(zi_i)M ■ ■ ■ . (11) 

The integrals are path ordered, 

•■■Zi+i > Zi > Zi-i---- (12) 

When the two operators A = 6/6E in He act on the two E in 
equation (0), the [E|]~^ must act on the holonomy between the 
E, because of the pattern of epsilons in He . Another pattern of 
epsilons would have forced thejE|]~^ to act on the holonomy to the 
right or to the left of the two E. 

I have chosen the pattern in equation (|I0|) because of its sym- 
metry and simplicity. The calculations of the next section are not 
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unduly sensitive to this choice: one can still find solutions when the 
[E|]^^ acts on the holonomy to the left or the right, although they 
differ slightly from the solutions when the [E^]"^ acts in the middle. 
Note that another symmetric choice yields nothing new: reversing 
the sign of the epsilons in equation (|^) does not change anything, 
since He is even under interchange of the A's. However, I have been 
unable to find a solution using an average over left, right, and middle 
holonomies, rather than the middle holonomy alone. 

I have verified that the pattern of epsilons given in equation ([l0|) 
is preserved by the constraint algebra, in the following sense. The He 
term is contained in the scalar constraint Hg. Consider a constraint 
commutator where Hg occurs on both sides of the equation. 

[Hs,H,]cxHg. (13) 

If the He on the left is given the pattern of epsilons shown in equa- 
tion (|TOD then the He on the right will also follow the same pattern. 
(Also, there will be additional terms of order epsilon on the right, 
because the point splitting violates diffeomorphism invariance; but 
these terms disappear in the limit e^O.) If there is no Hg on the 
right, 

[Hs,Hg]ocH,; 
[Hg,HG] = 0, 

then the right hand sides will contain additional terms of order ep- 
silon, which again disappear in the limit e— s-0. 

3 Solutions 

This section proposes a solution, then verifies that it is annihilated 
by the constraints. In II the constraints were expressed in terms of 
fields which are eigenstates of the surviving gauge invariance 0(2) 
or U(l) generated by Hg- That is, I used basis fields 

E| = [E^±iE^]/y2,a=x,y, (14) 

rather than the usual Ex and Ey; similarly, I used A^. The solutions 
of both II and the present paper are much simpler when expressed in 
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terms of these basis fields. Accordingly, I return to the Hamiltonian, 
equation and break it up into eigenstates of 0(2) by writing out 
the components of the Levi-Civita tensor, 

e-+ = -e+- = i, (15) 

while being careful to contract every + index with a - index, for 
example eMNA^B^ = e^+A+B" + ■ ■ -. 

Ht = ^(N'-N^)[-e,,,E-E°(E|)-ieedA+Aj-2iEb(9, + iA^)A+] 

+ 1(N' + N^)[-e^,E-E- (E|)-ieedA+A, + 2iEb (9, - iA,^)A^] 

-iNG[9.E| - i(E;A- - E^_A+)] + S.T. (16) 

It is perhaps also worth repeating the quantization equations using 
this 0(2) eigenbasis, since the pattern of ± signs may look unfamil- 
iar. 

At = h<5/5E- 

E| = -h6/6Al (17) 

To understand the pattern of signs, note that the two dimensional 
Kronecker delta in equation (^ has only off-diagonal elements when 
expressed in terms of 0(2) eigenstates: 6±^ = +1.) 

The ansatz for the solution contains n factors of E!^S+ followed 
by an equal number of factors of Ei!JiS_, where S± are SU(2) raising 
and lowering operators for the 2j + 1 dimensional representation of 
SU(2). The ES factors are separated by the z-axis holonomies M 
introduced at equation (||). 

^ /■Z2n+l 

ip{2n;i) = [[ dzi6'(zi+i - Zi)6'(zi,Zo) X 

i=l "'^0 

XM(z2n+l,Z2n)E;^'^(z2n)S-M(z2n,Z2n-l) " " ' E^"+' (z„+i) S„M(z,+i , Z,) X 

xE!l"(zn)S+M(z„Zn_i)E""-^(z„_i)S+---E!_^(zi)S+M(zi,zo) x 

xM(zo,Z2n+l)ea„+ia„- (18) 

If the holonomy M(zi+i, Zj) is pictured as a flux line extending from 
Zi to Zi+i, then ip corresponds to an open-ended flux hne stretching 
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from zq to Z2n+i, with 2n E operators inserted along the hne. The 
final M(zo,Z2n+i) holonomy represents a flux line which loops back 
from Z2n+i to Zq and turns the open-ended flux line into a closed 
loop. The SU(2) structure of equation (JT^) supports this interpre- 
tation. The S± and M are (2j+l) x (2j + l) dimensional matrices 
with rows and columns labeled by the eigenfunctions m^ of S^. The 
m^ subscripts on these matrices have been suppressed for simplicity, 
but clearly , since there are an equal number of raising and lowering 
operators in the chain, the flrst and last holonomies in the chain 
have the same value of m^. This means that the return segment 
M(zo, Z2n+i) has the right m^ subscripts to match the m^ of the flux 
exiting from the beginning and end of the open chain. Note that it is 
not necessary to sum over the mz (it is not necessary to take a trace) 
because the internal gauge group SU(2) has been flxed to 0(2), and 
the irreducible representations of 0(2) are all one-dimensional, la- 
beled by their m^ value. Later I use identities such as 

[Sz,S±] = ±S±, (19) 

but even here there is no need to sum over the intermediate values 
of m^ on the left, since only one intermendiate value occurs anyway. 
If I wished, I could replace every Si by the corresponding matrix el- 
ement; e.g. replace by m^; but that would be awkward, especially 
for the S±. 

Since m^ is not summed over, the m^ value at (say) the start of 
the chain remains a free parameter for the moment. Later in this 
section we shall see that the scalar constraint is not satisfied unless 
the holonomy M(zn+i, Zn) has m^ = —1/2. This particular holonomy 
occurs at the changeover point, where the chain shifts from S+ to 
S_. Note also the ai indices on either side of the changeover point 
must be antisymmetrized by the final epsilon tensor in equation (|18|). 
Again, this is required by the scalar constraint. The free parameters 
are therefore n, j, and the a^ away from the changeover point. 

Since every m^ occurs contracted with another m^, and every S± 
is contracted with an E^, it should be clear that ip is annihilated 
by the Gauss constraint. To verify this explicitly, note from equa- 
tion ( [TB| ) that this constraint is the sum of three terms. 
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The first term acts only on the holonomies: 

9,E|(z)M(zi+i, Zi) = iS,[5(zi+i - z) - 5(zi - z)]. (21) 

The remaining two terms act on the E^. These two terms, hke the 
first term, also generate factors of ±i5(zi — z), but no factors of S^. 
Now group together all terms in B.Qip containing a factor S{zi — z). 
If Zi is not_ an endpoint zq or Z2n+i, there will be three such terms, 
one from , and two from the holonomies on either side of E^ : 

i<5(zi - z)M(zi+i, Zi){±S^ + [S„ S^]}E^M(zi, z,_,) (22) 

From the commutation relations obeyed by the Si, the curly bracket 
vanishes. If Zi is an endpoint zq or Z2n+i, there will be no contri- 
butions from E fields, only contributions from the three holonomies 
in ip which depend on the endpoints. But these holonomies can be 
combined so that all dependence on the endpoints disappears: 

M(z2n+1, Z2n) " " " M(zi, Zo)M(zo, Z2n+l) = M(z2n, Zl)- (23) 

Hence there is no contribution from the endpoints either, and the 
Gauss constraint annihilates ip. 

Now consider the diffeomorphism constraint H^. The proof that 
this annihilates ip proceeds along the same lines as the corresponding 
proof for the solutions presented in II, since again i/; is a product 
of diffeomorphism- invariant integrals of the form / dzi (scalar den- 
sity E). I mention only the one point of the proof which is not 
straightforward. Since the range of integration is zq to Z2n+i, rather 
than the whole real axis, the points zq and Z2n+i are singled out. 
This seems to violate diffeomorphism invariance. However, the pre- 
cise statement of this invariance requires that the wavefunctional 
be annihilated by the smeared constraint / ^N^H^dz, where (5Nz is 
an infinitesimal change in the shift vector. Away from boundaries 
(5Nz is arbitrary, but at boundaries it must vanish, in order to pre- 
serve the boundary conditions N^ = 0. This vanishing is enough to 
guarantee the vanishing of any contributions to / 6NzB.zdzip from 
boundary points, and the proof of diffeomorphism invariance is not 
affected by the finite range of integration. 

Only the scalar constraint remains to be considered. From equa- 
tion ([l6| ) this can be written as a sum of three terms. 

Hs = H+ + H_ + He, (24) 
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where 

H± = ±iN'E^(9,TiA,")A^; (25) 

He = N'p)E(z)6ABeabA^(z + e)(E|(z))-iA^(z - e)/2. (26) 

It is understood that the and E| operators in these equations are 
to be replaced by the functional derivatives given at equations (|17D . 
Consider the action of }l± first. This operator acts on the factors of 
E^ in ip but does not remove the E; instead, it removes one of the 
adjacent holonomies. 

H±V = H±[- ■ ■ ^(zi+i - zO^(zi - Zi_i)M(zi+i, Zi)E5!(zi)S^M(zi, Zi_i) ■ ■ ■] 
= E[---^(2i+i-^i)^(2i-2i-i)M(zi+i,Zi)(±iN'E^(z)) X 

i 

X (9, T iA^)S{z - Zi)S^M(zi, Zi_i) ■ ■ ■]. (27) 

Now change the dz on the delta function to a —d^.^ and integrate 
by parts on Zi. The surface terms at z = zq, Z2n+i again vanish 
because the smearing function N' is actually a 6N' which vanishes at 
boundaries. The d^-^ acts on the M and the 6 factors in equation (^7h. 



Since The action on the M's brings down factors of A^, group these 
terms with the term in equation ( ^7] ) which already contains an A^: 

(5,,TiA,2)[M(zi+i,Zi)(±i)S^M(zi,Zi_i)] cx iA,^) [-iS,S^ + S^iS, ^ S^] 

= 0. (28) 

Thus the A^ term in equation (^) and equation (|27|) cancels terms 
where derivatives act on factors of M, leaving the terms where 
derivatives act on the 6 functions. 

H±^ = ±iN'^[---5(z-zO(^(zi+i-z05(zi-Zi_i)-5(zi+i-Zi)^(zi- 

i 

xM(zi+i,zOE^(z)S^M(zi,Zi_i)---]. 

The two delta functions in each term reduce one of the M's to unity, 
in effect removing it completely. As stated previously, the H± terms 
do not remove an E, but do remove a neighboring holonomy. 

Now consider the He term in the scalar constraint. From equa- 
tion (PBD, this removes two E and replaces them by an The two 
E which are deleted must be adjacent along the chain; otherwise 
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there will be two delta functions (5(z — Zi)5(z — zj) which will squeeze 
the range of integration of a dz^ to zero, where z^ lies between Zi 
and Zj. Because of the eAB in equation (^), one of the adjacent E 
must be an E^, and the other must be an E^. Therefore He acts 
only at the changeover point, where the chain shifts from S+ to S_. 

He^ = ■ ■ ■ Nf 2)E(z)e+_ea„+,a„5(z + e - z,,+{)5{z - e - z,)/2] x 

x(E|(z))-^S_M(z,+i,Zn)S+---ea,+,a„ 
= ■ ■ -Nf 2)E(z)e+_(e,.^^^,j25(z + e - z,+i)5(z - e - z^)/2] x 
xS_(-im,)-iM(z,+i,z,)S+--- 

I have used equation for the action of E^iz)'^ on M. Now take 
the limit e^O. 

He7/> ^ •■■N'p)E(z)(_i)(+2)5(z-z„+i)5(z--z„)/2] X 
xS_(-im,)-i(l)S+--- 
= • • • N'p)E(z)<5(z - z„+i)5(z - -Zn)]S_K)-iS+ ■ • {30) 

As predicted in the last section, the z±e point splitting is not needed 

to regulate ((5(z Zn))^ divergences; there are none. Rather, it is 

needed to define unambiguously which M the E|(z))~^ is to act on: 
it acts on the M between the two E grasped by the A's in He- 

The action of B.± and He on has now been described; it is 
time to put the pieces together and describe the action of the entire 
scalar constraint on ip. Each term in H^ip contains a double delta 
function (5(z — Zi+i)(5(z — Zi); cf. equation (^) and equation (|30D . It 
is therefore natural to group together terms having the same double 
delta function. Terms where Zj+i or Zi is a boundary point can 
be dropped, because N'(z) (or more precisely, 5N'(z) ) vanishes at 
boundary points. For the moment ignore the 6{z — Zn+i)5(z — Zn) 
terms at the crossover point, since these are the only terms to receive 
a contribution from the He piece of Hs. Hsip contains only two terms 
proportional to 6{z — Zi+i)5(z — Zj), i ^ n. One of the two comes 
from the term shown explicitly in equation (|29|). This term is the 
ith term in a sum over 2n terms. The other 6{z — Zi+i)5(z — Zi) term 
comes from the (i+l)st term in the same sum. The two terms have 
opposite sign and so cancel. 
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Now consider the 5(z — Zn+i)5(z — Zn) terms, those at the crossover 
point. Since the crossover point is flanked by both an E+^^ (at z = 
Zn+i)and an E!l" (at z = Zn), there will be contributions from both 
H_|_ and H_. For the H+ contribution, set i = n+1 in equation (P^; 
for the H_ contribution, set i = n. These two contributions have 
the same sign and do not cancel. 

(H+ + H_)7/; = ■■■N'5(z-z,+i)5(z-z,)2iE^^+^(z)S_E^_"(z)S+---ea„+,a„ 
= ■■■5(z-z„+i)5(z-Zn)2(2)E(z)S_S+---. (31) 

This will be canceled by the He contribution, equation (|30D , pro- 
vided one chooses 

m, = -1/2. (32) 

This must be the value of m^ at the changeover point. At the z = zq 
end of the chain, then, m^ = —1/2 — n. The n S+ operators raise 
this to a maximum of -1/2 at the changeover point; then the n S_ 
operators lower this back to m^ = — 1/2 — n at the z = Z2n+i end of 
the chain. At no point does m^ pass through the forbidden value of 
zero. This completes the proof that if), equation (0), is annihilated 
by all the constraints. 

From this solution one can generate others. E. g. by interchang- 
ing plus and minus SU(2) subscripts in equation (|18|) , one generates 
the "complex conjugate" wavefunctional, which is also a solution 
with all m^ positive and equation (^) replaced by m^ = -1-1/2. 
Also, it is possible to relax the requirement that all S± and M have 
the same j; see the concluding section of II. 

4 Discussion and Directions for Further Work 

What is the physical interpretation of the solutions of section 3? In 
II I argued that, since Hs generates time translations and Hz gener- 
ates space translations, it is reasonable to assume that the operators 
Hs±Hz represent displacements along the light cone directions t'±z' 
intrinsic to the plane wave metric. (In technical terms the metric 
considered here posseses two hypersurface orthogonal null vectors; 
and the corresponding hypersurfaces can be parameterized by t'±z' 
= constant, where the coordinates (t',z') are related to (t,z) by an 
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appropriate gauge transformation.) As a check, one can show that 
the commutator [Hs + Hz,Hs — H^] vanishes, as it should if it is 
the commutator of two independent translation generators. Note 
(Hs ± Hz)/2 is the operator multiplying N_' ± in equation (|TB|), 
and this operator is linear in = h5/(5Ej_. The solutions in II 
depend on either or fields, but not both. This implies that 
the solutions obey either (Hs + Hz)?/^ = or (Hs — Hz)?/^ = 0, which 
suggests that they depend on either t' — z' or t' + z', but not both. 
If the solutions of II represent radiation, then that radiation is uni- 
directional, traveling in a single direction along the z axis. The 
solutions constructed in the present paper contain both and E!^ 
fields. By the same argument, if the solutions of this paper repre- 
sent radiation, it is not unidirectional. These solutions correspond 
to scattering. 

The argument just outlined is very limited. It does reveal the di- 
rection of motion of the field. By itself, however, it cannot establish 
that the field is radiative, much less give details of the polarization 
and amplitude of the radiation. One can ask to what extent it is 
possible to go farther and construct operators which explore these 
details. In the 1970's, a great deal of work was done on charac- 
terization of classical radiative plane wave solutions, and most of 
this work is readily translated into the Ashtekar language. How- 
ever, it is one thing to characterize classical radiation, and another 
to translate a classical criterion into a quantum criterion. 

It is not hard to see what issues can arise. For simplicity, consider 
the corresponding problem in electrodynamics: given a radiation 
criterion valid for classical Maxwell electrodynamics, translate it 
into the corresponding criterion for QED. Consider, for example , 
the classical criterion that a plane wave traveling in the z direction 
has only a single polarization: 



where ip is now a state in the Hilbert space of the operators E and 
B. To obtain the quantum analog of the classical equation equa- 




(33) 



The corresponding quantum criterion is not 



(E, + iBy)ij = 0, 



(34) 
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tion (|33|), one must pass to large occupation numbers and use co- 
herent states: 

< coh I (Ex + iBy) I coh >= 0. (35) 

equation ( |5^ is far too strong: the vacuum state with respect to 
a given polarization need not be annihilated by the corresponding 
operator, equation (^5l), however, is 

too special, since it works for coherent states but not for eigen- 
states of the number operator. 

If one has the full machinery of QED, of course the correct proce- 
dure is to break up E and B into creation and annihilation operators 
and demand that aip = 0, a. an annihilation operator for the relevant 
polarization. In a generally covariant theory, creation and annihi- 
lation operators are not automatically available, since normally a 
preferred time coordinate is required for their definition, e. g. the 
time associated with the timelike Killing vector in a static space- 
time. Here one has the preferred coordinates t' ± z', but in these 
coordinates the metric is conformally flat rather than static. On 
the other hand, unidirectional plane waves are known to obey a su- 
perposition principle |TT|, |T^, which suggests that the unidirectional 
case may possess particle-like modes. Even if no quantum criteria 
are forthcoming, the classical criteria by themselves are quite useful 
and illuminating, and I plan a paper which will (at a minimum) 
translate the classical criteria for plane radiation into the Ashtekar 
language. 
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